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ABSTRACT 

We find statistically significant correlations in the cosmological matter power spec- 
trum over the full range of observable scales. While the correlations between indi- 
vidual modes are weak, the band-averaged power spectrum shows strong non-trivial 
correlations. The correlations are significant when the modes in either one or both 
bands are in the non-linear regime, and approach 100% for pairs of bands in which 
all the modes are non-linear. The correlations are weaker, but not absent, when com- 
puted in redshift space. Since estimates of the power spectrum from galaxy surveys 
require band-averaging, the correlations must be taken into account when comparing 
a measured power spectrum with theoretical models. 
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1 INTRODUCTION 

In our current paradigm for cosmological structure forma- 
tion, initially small, primordial density perturbations grow 
through gravitational instability to form the large-scale 
structures that we see in the distribution of galaxies today. 
The most popular assumption is that the primordial fluctua- 
tions were distributed according to a homogeneous Gaussian 
random process. Under this assumption all of the statisti- 
cal information is encoded in the power spectrum. The as- 
sumption of initially Gaussian modes is a prediction of the 
simplest models of inflation, and has received some limited 



(Kogut ct al. 1996; 


Heavens 1990 and large-scale structure 


(Bouchct ct al. 1993 


: Gaztanaga 1994; Nusser, Dekel & Yahil 


1994 Feldman, Kaiser & Peacock 1994; Colley 1997J) . How- 



ever, it is known that as the perturbations grow and become 
non-linear, the modes become coupled. In this paper, we 
study this process quantitatively using iV-body simulations. 

Since it is the lowest order non-vanishing statistical de- 
scription of the density or galaxy field, the power spectrum 
has been at the focus of much recent attention. Many galaxy 
surveys have attempted to determine the nature of the power 
spectrum in two and three dimensions, and two massive sur- 
veys, the Anglo- Australian Two Degree Field (2dFQ) and 

the Sloan Digital Sky Survey (SDSS^) are moving the sub- 
ject toward high-precision measurements. Several methods 



* http:/ /meteor. anu.edu.au/~colless/2dF 
t http:/ /www. astro. princeton.edu/BBOOK 



for performing detailed comparisons between a measured 
power spectrum from a galaxy survey and a theoretical 
power spectrum have been devised (see Tegmark et al. 1998 
for a review). Critical to assessing the accuracy of the mea- 
surements is an evaluation of the statistical properties of the 
power spectrum. To date most estimates of the precision to 
which the power spectrum may be measured have assumed 
that the Fourier components of the density fluctuations are 
still distributed as a Gaussian random process. In that case 
the standard deviation of the power spectrum is determined 
by the value of the power spectrum. We expect this assump- 
tion to be valid at early times and on sufficiently large scales 
for models with Gaussian initial conditions. The question of 
how large is "sufficiently large" requires a detailed calcula- 
tion. 

We show here that non-linear clustering gives rise to sig- 
nificant correlations in the band-averaged power spectrum 
at different scales in the regime of observational interest. 
The correlations increase with increasing (spatial) frequency, 
reaching levels near 100% when the perturbations are non- 
linear. Moreover, the correlations require the modes in only 
a single band to be non-linear: we find significant correla- 
tions between high and low frequencies even if the modes 
in the low frequency band are still in the linear regime. A 
consequence of the mode coupling is an increase in the dis- 
persion of the power spectrum, suggesting that estimates of 
the accuracy to which the power spectrum may be measured 
should be revisited. 

The outline of the paper is as follows: in §^ we discuss 
the expectations for correlations in the power spectrum be- 
tween modes. In §^ we describe our numerical experiments, 
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while in §^ we discuss the statistical properties of the corre- 
lations. We finish in §^ with some comments on the impact 
of our results for future surveys designed to measure the 
matter power spectrum. 



2 CORRELATIONS IN THE POWER 
SPECTRUM 

The density contrast is defined by <5(r) = [p(r) — p}/p, where 
p(r) is the matter density at position r, and p is the mean 
density. We consider a large spatial volume V and impose 
periodic boundary conditions, thus r is continuous but its 
Fourier conjugate k is quantized. One can take the con- 
tinuum limit by replacing V~ 1 ^ k with J d 3 k/(2ir) 3 . Our 
Fourier transform convention has 



5k = V~ 



and 



5{v) 



V 



(2tt) 3 



d r 5(v) exp[ik • r] 



d k <5k exp[— ik ■ r] 



(1) 



(2) 



both dimensionless. We define the power spectrum as 
P(k) = V (|<5fc| 2 ) — 1/fl, where the angled brackets {• • •) de- 
note an ensemble (and Poisson) average, and the shot noise 
term 1/n has been subtracted. Here, n is the mean num- 
ber density of objects used to measure the power spectrum. 
The power spectrum is related to the 2-pt spatial correlation 
function £(r) by 



P(Jfc) = / d 3 r £(r)exp[ik-r] 



(3) 



We denote our (unbiased) estimate of the power spectrum 
at k, which we obtain from the simulations, as 



P(k) = V\S k \ 2 -1/n 



(4) 



Typically we reduce the scatter in this quantity by averaging 
over a thin shell in k-space with |k| ~ k. 

The covariance between two random variables x and 
y is defined as cav[x,y] = (xy) — (x)(y). The variance is 
var(:r) = cov[x,x]. The correlation between x and y is de- 
fined as p(x,y) = cov[x,y]/[vai(x)vaT(y)] 1 ^ 2 . The covari- 
ance of our power spectrum estimator, P(k), depends on 
the 2-pt function (power spectrum), the 3-point function 
(bi-spectrum) and the 4-point function (tri-spectrum). We 
obtain, to leading order in (nV) -1 as nV — > oo (Peebles 
1980, §36), 



cov[P(k),P(k')] = [P^+n- 1 ] (5 k ,k' + 5 k ,_ k 

+ r(k,-k J k' 1 -k')/v 



(5) 



with <5k,k' = 1 for k = k', and otherwise, and where we 
have written the tri-spectrum as T. The tri-spectrum is de- 
fined as the Fourier transform of the 4-point function n ac- 
cording to 



T(k 1 ,k 2 ,k 3 ,k 4 ) = 



— / d 3 xi d i X2 d A X3 d 3 X4. r?(xi, X2, X3, X4) x 
exp[iki ■ xi + ik.2 ■ X2 + ik 3 ■ X3 + ik.4 ■ X4]. 



(6) 



The first term in Eq. (|^) is the usual result for Gaussian fluc- 
tuations, viz var[a; 2 ] = 2(x 2 ) 2 for a real variable x. Because 
<5 k is complex with uncorrelated real and imaginary parts, 
P = 2([Re(<5 k )] 2 ) = 2([Im((S k )] 2 ), and the Gaussian contri- 
bution to the error on P becomes equivalent to the total 
power spectrum: the sum of the signal and the noise. Only 
the tri-spectrum contributes to the off-diagonal (k 7^ k') 
elements of the covariance. 

The terms sub-dominant in nV, all shot noise terms, 
are also straightforward to derive. There is a constant term 
(Vn 3 ) -1 , contributions from the 2-pt function 

-pL [P(]k + k'|) + P(|k - k'|) + 2P(k) + 2P(k')] (7) 
and the 3-pt function 

' [P(k,-k,0) +B(0,k',-k') + 

P(k + k', k, -k') + P(k - k', -k, k') + 
P(k,k'-k, -k')+P(k, -k-k',k')] (8) 



hV 



where the bi-spectrum is defined as the Fourier transform of 
the 3-pt function £ 



B(ki,k 2 ,k 3 ) = — J d 3 X! d s x 2 d 3 x 3 C(xi,x 2 ,x 3 ) x 

exp [iki ■ xi + ik 2 ■ x 2 + ik 3 ■ x 3 ] . (9) 

Finally we consider the result for band-averaged esti- 
mates of P{k). Denote by Pi our estimator P(k) averaged 
over a set of Ni k;, a with |k; ja | ~ fc; for a = 1 to Ni. 
Straightforwardly 



1 



NiN„ 



EE 

a = l /3 = 1 



P(k l , a ),P(k 



3,0) 



.(10) 



Even when the bands i and j are disjoint, we shall see that 
the band-averaging introduces non-trivial correlations in the 
power spectrum. These correlations are given by the tri- 
spectrum averaged over the configurations in the shell. In 
principle one could apply different weights to the modes in 
the band, but we have not pursued this line of inquiry. It is 
traditional to merely average over the directions in a fc-shell 
and this is what we have chosen to do in the simulations. 



3 THE GROWTH OF CORRELATIONS 

Gravitational perturbation theory (PT) has given some in- 
sight into the growth of clustering beyond the linear regime, 
and in particular the effect of clustering on the evolution of 
t he power spe c trum, the bi-spec t rum, and the t r i-spectrum 
( |Peebles 198o| ; |juszkiewicz 198l| ; |Vishniac 19"83j ; [Fry 1984j ). 
Perturbation theory shows that the variance in the power 
spectrum will grow for modes that are still well in the linear 
regime. In addition, PT shows that the non-linear gro wth 
of de nsity fluctuations will give rise to a tri-spectrum ( fcry| 
1984). The form the tri-spectrum takes in PT involves mul- 
tiples of the first order power spectrum and no powers of 
V. While this is shown explicitly by Fry to the lowest non- 
vanishing order in the tri-spectrum, it is straightfoward to 
extend the result to all orders. Consequently, the contribu- 
tion of the tri-spectrum to the power spectrum covariance 
between any off-diagonal pair of k- modes is of order 1/V in 
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PT, and so vanishes in the lar ge vol ume limit according to 
equation (§). Fan & Bardeen ( |l995| ) have argued that the 
power covariance matrix should remain diagonal even in the 
non- linear regime. 

Nonetheless, as we describe below, correlations in the 
band-averaged power spectrum may remain even in the limit 
V — ► oo. This is a consequence of the decrease in the vari- 
ance of the power spectrum with the increased number of 
modes in the band. Since the number of modes in a band 
is proportional to V, the 1/V suppression of the covariance 
between modes cancels in the expression for the correlation 
p, leading to a finite value even as V — * oo. We investigate 
the development of the correlations in the non-linear regime 
using iV-body simulations. 

We remark that for a finite-volume survey, we would 
expect the presence of both covariances and correlations be- 
tween individual modes according to equation (|s|) and the 
associated shot-noise terms. 

Using a particle-mesh (PM) code, described in detail 
in Meiksin, White & Peacock (1999), we have performed 



a series of several thousand realizations of the evolution of 
clustering in a ACDM universe. The parameters chosen are 
fto = 0.4, n A = 0.6, h = 0.65, n B h 2 = 0.03, and n = 1.030, 
where n is the primordial spectral index. The power spec- 
trum is COB E-normalized using the method of Bunn & 
White (1997). The slight tilt has been chosen t o reproduce 



the ab undance of rich clusters of galaxies today ( White, Efs- 



tathiou fc Frenk 1995 ), specifically the rms mass fluctuation 
in an 8 h~ Mpc sphere is erg ~ 0.92. The non-linear scale, 
where the variance becomes unity, is fc n i = 0.19 h Mpc -1 . 
Tests involving the higher order moments indicate that the 
code should accurately r eproduce th e bi- and tri-spectrum 
on the scales of interest (White 1998). 



All the simulations were run with a 128 force grid 
and either 128 3 or 256 3 particles to isolate the effects of 
shot noise. Two box sizes were used, 200 /i -1 Mpc and 
400 h^ 1 Mpc , to test for finite volume effects. For each box 
size, N ~ 10 3 realizations with different Gaussian initial 
conditions were evolved from 1 + z = 20 to the present. For 
the 256 3 simulations, N ~ 10 2 realizations were performed. 
We verified that starting the simulations at 1 + z = 30 did 
not alter our results. For each realization the power spec- 
trum was calculated in 20 bins logarithmically spaced in k 
from the fundamental mode to the Nyquist frequency of the 
force grid, as described in Meiksin et al. (1999). The en- 
semble averages of the previous section were approximated 
by an average over the N realizations. Convergence of the 
correlations was typically obtained with several hundred re- 
alizations. 



The resulting correlations in the real space power spec- 
trum at z = for several selected frequency bands from 
the 200 /i -1 Mpc simulations are shown in Table [j]. For 
each entry we may assess the likelihood of no correlation 
using the Spearman rank order correlation coefficient r 3 
( |5iegel fc Castellan 1985] ). The probabilities are < 0.001 
for pairs of bands with both central frequencies exceeding 
0.14 ftMpc -1 . An estimate of the probability for all the 
correla tions is provided by the Kendall concordance statis- 
tic W (siegel & Castellan 1985), which is linearly related to 



0.6 



0.4 



0.2 



- 



Correlations 
(200 h" 1 Mpc box) 



I I I 



_L 



0.1 

k (h Mpc -1 ) 



Figure 1. The correlations in power between a band centred at 
k and bands centred at k = 0.031/iMpc (triangles) and k = 
0.17/i Mpc - 1 (squares). The correlations are shown both in real 
space (filled) and redshift space (open). The correlations become 
significant when the modes in one of the bands enter the non- 
linear regime. Redshift space distortions reduce the correlations 
on small scales. 



the average of the r B values for all the band pairs^j For the 
12 frequency bands shown in Table [l], we obtain W = 0.33, 
with a vanishingly small probability for obtaining a value 
so large assuming no correlations. Indeed, it is only for cor- 
relations for which both bands are confined to frequencies 
k < 0.093 h Mpc -1 that the collective probability for non- 
correlation is found to exceed 0.001. 

We also compute the correlations in redshift space by 
including the peculiar velocities of the particles. The cor- 
relations are weaker (W = 0.26 for the 12 bands), but are 
still highly significant. We expect that the correlations are 
weaker in redshift space both because of the reduced power 
on small scales in redshift space and because the peculiar 
velocities introduce extra randomness which destroys the 
correlations. 

In Figure [j], the correlations are shown between 
each band centred at k and the bands centred at k = 
0.031 ft Mpc -1 or k = 0.17 /iMpc -1 . The error bars indi- 
cate one standard deviation calculated from a p — N^ 1 ^ 2 , 
assuming the modes are drawn from a Gaussian distribu- 
tion.^] In the last row of Table [l], we show the density fluc- 
tuations on each scale k, as defined by a = [A 2 (k)] 1 ^ 2 , 

t The relation between W and the average coefficient f s for Nf, 
bands is W = (1 — N~ 1 )f s + N~ 1 . The correlation coefficient r s 
between two quantities is the product-moment coefficient com- 
puted by replacing the values of the quantities by their ranks. 
S The error is independent of the number of modes entering the 
band averages. The errors will be reduced in the presence of corre- 
lations. The error bars shown demonstrate the significance of the 
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Table 1. Power spectrum correlation matrix for 200 h 1 Mpc simulation. The last row shows cr = [A 2 (fc)] 1 / 2 . 



k (TiMpc- 1 ) 


0.031 


0.044 


0.058 


0.074 


0.093 


0.110 


0.138 


0.169 


0.206 


0.254 


0.313 


0.385 


0.031 


1.000 


-0.017 


0.023 


0.024 


0.042 


0.154 


0.176 


0.188 


0.224 


0.264 


0.265 


0.270 


0.044 


-0.017 


1.000 


0.001 


0.024 


0.056 


0.076 


0.118 


0.180 


0.165 


0.228 


0.234 


0.227 


0.058 


0.023 


0.001 


1.000 


0.041 


0.027 


0.086 


0.149 


0.138 


0.177 


0.206 


0.202 


0.205 


0.074 


0.024 


0.024 


0.041 


1.000 


0.079 


0.094 


0.202 


0.229 


0.322 


0.343 


0.374 


0.391 


0.093 


0.042 


0.056 


0.027 


0.079 


1.000 


0.028 


0.085 


0.177 


0.193 


0.261 


0.259 


0.262 


0.110 


0.154 


0.076 


0.086 


0.094 


0.028 


1.000 


0.205 


0.251 


0.314 


0.355 


0.397 


0.374 


0.138 


0.176 


0.118 


0.149 


0.202 


0.085 


0.205 


1.000 


0.281 


0.396 


0.488 


0.506 


0.508 


0.169 


0.188 


0.180 


0.138 


0.229 


0.177 


0.251 


0.281 


1.000 


0.484 


0.606 


0.618 


0.633 


0.206 


0.224 


0.165 


0.177 


0.322 


0.193 


0.314 


0.396 


0.484 


1.000 


0.654 


0.720 


0.733 


0.254 


0.264 


0.228 


0.206 


0.343 


0.261 


0.355 


0.488 


0.606 


0.654 


1.000 


0.816 


0.835 


0.313 


0.265 


0.234 


0.202 


0.374 


0.259 


0.397 


0.506 


0.618 


0.720 


0.816 


1.000 


0.902 


0.385 


0.270 


0.227 


0.205 


0.391 


0.262 


0.374 


0.508 


0.633 


0.733 


0.835 


0.902 


1.000 


a 


0.17 


0.25 


0.33 


0.44 


0.55 


0.62 


0.76 


0.91 


1.07 


1.29 


1.55 


1.87 



where A 2 (fc) = k 3 P(k)/(2n 2 ). While the correlations be- 
tween pairs of low frequency bands, both well in the lin- 
ear regime (k < fc n i), are consistent with zero, the power 
spectrum correlation between low frequency bands and high 
frequency bands rapidly rises as the frequency of the latter 
increases. The correlations involving the higher frequency 
bands are significantly reduced in redshift space. 

In Figure ^, we show that no significant difference is 
found between the 200 /i -1 Mpc and 400 h" 1 Mpc box simu- 
lations, suggesting that the correlations have converged and 
are not due to the finite volume of the simulations. Similarly, 
we show in Figure ^| that the correlations do not decrease 
as the number of particles increases, demonstrating that the 
correlations are not dominated by shot noise contributions. 

As an additional test, we performed a similar set of 
~ 10 3 runs using the Zel'dovich approximation with the 
same box sizes and particle numbers as for the PM code. 
With the Zel'dovich calculations we could displace the par- 
ticles from random positions within the simulation volume, 
from a regular grid or from random positions near grid zones. 
This allowed us to test that our initial conditions did not 
contribute spuriously to the correlations. We found, see Fig- 
ure ^, that the correlations appeared much more strongly in 
the PM runs than the Zel'dovich runs, suggesting that they 
are induced by gravitational instability, as shown in Fig- 
ure^, and are not an artifact of our numerical technique. We 
note, however, that evolving the particles in the Zel'dovich 
approximation, even when the particles were initially placed 
completely at random within the simulation box, again re- 
sulted in statistically significant correlations. 

We have also tested the random number generation by 
using two separate generators for 512 Zel'dovich runs each. 



one ba sed o n a multiplicative congruential generator ( Pres; 
et al. '. 992 ), and the other based on a lagged Fibonacci 



generator ( Bertschinger 1995), the periods of both of which 
should be ample for our purposes. The results were statisti- 
cally indistinguishable. 



correlations, but do not strictly show the error in the calculated 
correlation strength when the correlations arc non-vanishing. 
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solid 400 h Mpc box 
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0.1 

k (h Mpc -1 ) 



Figure 2. The correlations in real space power between the bands 
centred at k and k = 0.17/iMpc — 1 (open circles) in the 200 
h~ 1 Mpc simulation. Also shown are the correlations for the k = 
0.16/1 Mpc -1 band (filled triangles) and the k = 0.19/iMpc -1 
band (filled circles) in the 400 Mpc simulation. The filled 
points are slightly offset in frequency for clarity. No significant 
reduction in the correlations is found in the larger box compared 
with the smaller, suggesting that the correlations are not a finite 
volume effect. 



4 THE STATISTICAL PROPERTIES OF THE 
CORRELATIONS 

To trace the origin of the correlations, we create two new 
bands with only 32 modes each for the 200 /i _1 Mpcbox 
simulation. One band is centred at k — 0.47 ft Mpc -1 and 
the other at k = 0.71 h Mpc -1 . We have first verified that 
the density fluctuations 8k themselves are uncorrelated, as 
they must be by homogeneity. We find W = 0.008, with a 
probability of non-correlation of p = 0.44. For the two band- 
averaged P(k) the correlation is p = 0.23, with a probability 
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Figure 3. The correlations in power between the band centred 
at k and the bands centred at k = 0.031ft Mpc -1 (triangles) 
and k = 0.17ft Mpc -1 (squares). The correlations are shown for 
the 200 ft, -1 Mpc box simulations for 128 3 (open) or 256 3 (filled) 
particles. 



0.8 -□ Zel k = 0.06 
- . pm k = 0.06 

0.6 Zel k = 0.17 
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Figure 4. The correlations in power between the band cen- 
tred at k and the bands centred at k = 0.06h Mpc -1 (squares) 
and k = 0.17ft Mpc -1 (triangles) computed using either the 
Zel'dovich approximation (open) or the PM code (filled). While 
the scatter at low-fe is large due to the small number of modes 
present, the PM correlations are stronger than the Zel'dovich cor- 
relations. 



p = 6.3 x 10 that the two are uncorrelated. For the full 
set of (64 x 63)/2 P{k) correlations, we find W = 0.025, 
with a probability p — 2.3 x 10 -22 , demonstrating the clear 
presence of highly significant correlations among the indi- 
vidual P(k). However, we find that none of the values in the 
32 x 32 inter-band correlation matrix for the individual P(k) 
estimates is as great as the correlation value p = 0.23 found 
between the band-averaged power spectra above. Indeed, the 
average correlation value in the matrix is 0.0099, and the 
standard deviation 0.037. We conclude that the strong cor- 
relations among the band-averaged power spectra are built up 
in the band-averaging procedure itself from the much weaker 
individual mode correlations. 

We may understand the origin of the correlations by 
considering a thin shell of width Afc centred at ki. Neglect- 
ing the shot-noise terms, the variance in the band-averaged 
power Pi(k) is var[Pj(fc)] = 2Pi(k) 2 /N k +T/V, where T de- 
notes an average over the configurations in the k— shell, and 
iVfc = 1/(27r) _3 47rfc 2 Ak is the number of modes in the shell. 
When the contribution of the tri-spectrum to the variance is 
negligible, we find that the correlation in the band-averaged 
power is given by p[Pi(k), Pj(k')] ~ [T /P(k)P 3 (k')]kk' Ak, 
and V has explicitly cancelled. Here T denotes the tri- 
spectrum averaged over all the configurations between both 
k— shells. For logarithmically spaced frequencies this be- 
comes p ~ [T/Pi(k)Pj(k')]{kk') 3/2 A\ogk. 

To make further progress, let us assume the hierarchi- 
cal clustering ansatz. The form of the reduced 4-point func- 
tion rj is then T](xi, x 2 , x 3 , x 4 ) = i? a [£(xi2)£(x23)£( x 34) + 
eye. (12terms)] + i?f,[£(xi2)£(xi 3 )£(xi4) + eye. (4terms)], 
where Xj j = Xj — x^ and R a an d Rt are constants of or- 
der one (Peebles 1980; Fry 1984). Using equation (n), we 
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Figure 5. The correlations in real space power between the bands 
centred at k and k = 0.17ft Mpc -1 . The correlations are shown 
for 3 different redshifts: z = 0.5 (circles), z = 0.3 (triangles), and 
z = (squares). The correlations grow with time. 



obtain for the contribution to the off-diagonal power spec- 
trum covariance 



V -1 T(k,-k,k',-k') 



(11) 



R 



^ [p(|k - k'|) + p(|k + k'|)] [p(k) + p(k')] + 
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Rb 
V 



P(k)P(k') [P(k)+P(k')] 



We find then that the power correlations for individual 
modes of frequencies k and k' will be on the order of 
[P(k) + P(k')]/V ~ 1(T 3 - 1(T 2 for a 200/1" 1 Mpc box. 
After band-averaging, and for P(k') ^> P(k), the cor- 
relations become p ~ Pj(k')kk' Afc for linearly spaced 
fc-shells, and p ~ A 2 (k')(k/k') 3/2 A log k for logarithmic 
spacing. Thus for logarithmically spaced shells, at a fixed 
k! the correlations will increase with k like k 3 ^ 2 . When 
A 2 (fc) ~ (Alogfc) -1 , the tri-spectrum term will domi- 
nate the variance of Pi and the correlations will flatten to 
p~ [A 2 (k)]' 1/2 A 2 {k')(k/k') 3/2 (Alogk) 1/2 . This is just the 
behaviour displayed in Figure hi Motivated by the hierar- 
chical ansatz form, we find the correlations for k S> k' and 
P(k') ^> P(k) are well approximated by 



(Ra + R 



A 2 (/fc')Alogfc 



R b [A 2 0) A log k 



1/2 : 



(12) 



where R a = (34/21) and Rb = 682/189 are angle- ave raged 



coefficients derived in perturbation theory ( Fry 1984 ). We 
find that a direct computation using equation (|ll|) with R a 
and Rb agrees less well with the simulations, over-estimating 
the correlations by as much as a factor of 2. It is unclear 
why the direct computation fails, but it improves toward 
the higher k shells, which may indicate that using the angle- 
averaged coefficients is an inadequate approximation for the 
bins with fewer modes. Since performing the full set of inte- 
grals in perturbation theory for all the modes in the k— shells 
is extremely cumbersome, we will not pursue the hierarchi- 
cal ansatz any further, but note that it may be worth greater 
consideration. 

A consequence of the correlations is a reduction in the 
number of degrees of freedom describing the variance of the 
power spectrum. While the distribution of band-averaged 
power is still x 2 , the modes are no longer independent and 
the "effective number" of degrees of freedom will be less 
than the number of modes in the band. We may define the 
"effective number" of degrees of freedom ridof by 



ridof 



2(P} 2 
var(P) ' 



(13) 



Since the density fluctuation 8k for each mode consists 
of independent real and imaginary parts, when the modes 
are Gaussian-distributed each mode will contribute two de- 
grees of freedom to the average. The band-averaged power 
spectrum will then be distributed like \ 2 with ridof degrees of 
freedom. Correlations between modes will reduce ridof below 
the number expected from mode counting. This is a strong 
effect at high frequencies, as shown in Figure []. Motivated 
by the hierarchical ansatz, we find that the reduction factor 
is approximately described by 1 + 2Rb A 2 (k) A log k. While 
fidof agrees with the number of degrees of freedom expected 
for Gaussian-distributed modes in the linear to quasi-linear 
regimes (k < 0.1 ftMpc -1 ), on scales for which the fluc- 
tuations are non-linear ridof falls substantially short. The 
errors on the power spectrum in the non-linear regime will 
correspondingly greatly exceed estimates based on the as- 
sumption of Gaussian statistics. 

In spite of the presence of strong correlations, the dis- 
tribution of band-averaged power spectra is still reasonably 
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Figure 6. The number of degrees-of-freedom expected in each 
frequency band based on counting modes in the simulation box 
(squares) vs the equivalent number ridof found in the simula- 
tions (circles). Also shown are maximum-likelihood estimates of 
n dof given by fitting the power spectra to a x 2 distribution 
(crosses). Correlations between modes reduce the effective num- 
ber of degrees-of-freedom at high frequencies, and so increase the 
dispersion in the band-averaged power spectrum estimates over 
the expectation for uncorrelated modes. 



well-described by a x 2 distribution. We show the maximum- 
likelihood estimates for the numbers of degrees of freedom 
in Fig. |^ found by assuming the power spectra in a given 
band are distributed like x 2 ■ Formally, the KS test rejects 
the x 2 distribution with 90-99.7% confidence at frequencies 
in the range 0.25 < k < 2 /iMpc . (We haven't examined 
higher frequencies.) The distribution of the redshift-space 
estimated power spectra, however, is found to be very close 
to x 2 at a h frequencies. 



5 CONCLUSIONS 

Estimates of the matter power spectrum from galaxy sur- 
veys have generally assumed the power in separate modes 
to be uncorrelated. The gravitational growth of the modes, 
however, will give rise to a tri-spectrum which will induce 
correlations in the band-averaged power spectrum. In a set 
of numerical experiments, we find that significant correla- 
tions in power develop when at least one of the modes is 
non-linear, and approaching 100% when the modes in both 
bands are non-linear. Although the power spectrum correla- 
tions are weaker in redshift space, they are still statistically 
significant. 

In addition to redshift space distortions, the measured 
galaxy power spectrum will also depend on the bias be- 
tween the galaxy count fluctuations and the underlying mat- 
ter fluctuations. We have not computed the correlations 
in power allowing for bias. Although this is straightfor- 
ward within any given bias scheme, which scheme, if any, 
the true galaxy distribution follows is unclear at this time. 
More direct measurements of the dark matter fluctuations 



may be obtained throug h galaxy velo city flows (Dckel 1994 
Strauss fc Willick 1995[) or lensing flBlandford et al. 1991 
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Miralda-Escudc 199lj kaiser 1992|; |Kaiser 19981 kcliak 1998| 



Bernardcau, van Waebeke fc Mellier 1997; Jain, Seljak & 



White '. 998 ) . Our results are not directly applicable to these 
measures, as they involve projections (integrals) of the mat- 



ter 



Hfc" 



po wei spectrum. However, we expect mat measme - me. 
errors will again exceed thooo cotimatod aaauming Gauaaian — 
distributed density fluctuations as a result of correlations in 
the matter power spectrum. 

The presence of power spectrum correlations restricts 
the range of applicability of many techniques for fitting 
models to the measured power spectra (e.g. Tegmark et al. 
1998) to very large scales. In particular, it may no longer 
be assumed that the fe-space density fluctuations 5k are dis- 
tributed as a Gaussian random process: the tri-spectrum 
must explicitly be taken into account if any of the modes fit 
are in the non-linear regime. This constraint may be partic- 
ularly severe if the window function of the survey is broad, 
mixing modes across a wide range of scales, and/or the sur- 
vey is not very deep. Correlations will affect the model- 
fitting based on all completed galaxy surveys, and may only 
be circumvented with wide-angle deep surveys like the Sloan 
Digital Sky Survey and the 2dF survey now underway. For 
a catalog the depth of the SDSS Northern Polar Cap red- 
shift survey, the minimum freque ncy spacing for unc orre- 
lated bins is Ak w 0.015 ftMpc" 1 flylciksin et al. 1999j ). On 
a scale of k = 0.19 h Mpc -1 , the scale of non-linearity for the 
model investigated here, the expected correlation with the 
power at k! = 0.1 ftMpc" 1 will be on the order of p ~ 0.15. 
Fitting a given model to measurements of the power spec- 
trum in the non-linear regime will require computing the cor- 
relations for the model. Provided there is a sufficient number 
of modes per band, the Central Limit Theorem guarantees 
the band-averaged power spectrum estimates will be dis- 
tributed as a multivariate Gaussian random process, simpli- 
fying the fitting process. Designing an efficient fitting pro- 
cedure in the presence of correlations is a topic worthy of 
future investigation. 
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